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Abstract Superfluid 3 He- A shares the properties of spin nematic and chiral or- 
bital ferromagnet. Its order parameter is characterized by two vectors d and 1. This 
doubly anisotropic superfluid, when it is confined in aerogel, represents the most 
interesting example of a system with continuous symmetry in the presence of ran- 
dom anisotropy disorder. We discuss the Larkin-Imry-Ma state, which is charac- 
terized by the short-range orientational order of the vector 1, while the long-range 
orientational order is destroyed by the collective action of the randomly oriented 
aerogel strings. On the other hand, sufficiently large regular anisotropy produced 
either by the deformation of the aerogel or by applied superflow suppresses the 
Larkin-Imry-Ma effect leading to the uniform orientation of 1. The interplay of 
regular and random anisotropy allows us to study many different effects. 
PACS numbers: 61.30.-v, 67.57.-z, 75.10.Nr. 



1 Introduction 

NMR experiments on liquid 3 He confined in aerogel demonstrate a rich but still 
unexplained life in two regions of the phase diagram in Fig. [TJ (i) in the region 
between the lines T c \, and 7" ca of superfluid phase transition in bulk liquid and in 
aerogel correspondingly 1 (in principle, Berezinskii odd-frequency pairing state 
is possible there 2 ); and (ii) in the region of the so-called A-like phase which is 
metastable in the absence of external magnetic field 3,4,5 . Here we concentrate on 
the A-like phase, which will be discussed in terms of the Larkin-Imry-Ma effect 
occurring in the anisotropic 3 He-A superfluid when it is confined in aerogel. 
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Fig. 1 Two main problems for superfluid 3 He in aerogel: (1) What is the state of 3 He between 
two lines of superfluid phase transition? What is the state of 3 He in the so-called A-like phase? 

1.1 Doubly anisotropic superfluid 3 He-A 

Superfluid 3 He-A shares the properties of spin nematic and chiral orbital ferro- 
magnet. Its order parameter is characterized by two vectors d and l 6 : 

A ili = Ae'^d^riii + ihi) , i = mxn. (1) 

Unit vector d marks direction of anisotropy axis in spin space, and is responsible 
for anisotropy of spin susceptibility; unit vector 1 marks the direction of the orbital 
angular momentum of Cooper pairs and simultaneously the direction of axis of the 
orbital anisotropy and it is responsible for anisotropy of superfluid density: 

F, + F 0+ F„ = ^ + ^ + M(y f .V + ^(y,xV- to (a.V. (2) 
2Z|| %X± 2 2 

Here S is spin density in the applied magnetic field, S a = Zap Hp ; v, is the velocity 

of superfluid mass current f s = p' s 3 v s . Because of the chiral nature of the orbital 
order parameter m + m, the superfluid velocity and its vorticity 

v s =^-(V0+m,Vn ( ) , Vxv s = ^e ij %VljxVl k , (3) 
2m 4m 

can be generated by the texture of the 1-vector 7 . The last term in Eq.© describes 
a tiny spin-orbit coupling F so between d and 1, which is however very important 
for the NMR measurements since it produces the frequency shift. 
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short-range order 



Fig. 2 While the short-range order is well defined, the collective orientational effect of aerogel 
strings destroys the long-range order on a macroscopic Larkin-Imry-Ma scale Llim 



This doubly anisotropic superfluid, when it is confined in aerogel, represents 
the most interesting example of a system with continuous symmetry in the pres- 
ence of random anisotropy. For such systems the surprizing conjecture was made 
by Imry and Ma 8 that even a weak disorder may destroy the long-range orien- 
tational order (LROO). This is the counterpart of the effect of collective pinning 
in superconductors predicted by Larkin 9 , in which weak impurities destroy the 
long-range translational order of the Abrikosov vortex lattice. Recent NMR ex- 
periments allow us to discuss the state of 3 He-A in aerogel in terms of the Larkin- 
Imry-Ma (LIM) state. The LIM state in 3 He-A is characterized by short-range 
orientational order of the orbital vector 1, while the LROO is destroyed by the 
collective action of the randomly oriented aerogel strings (Fig. 0. The order of 
magnitude theoretical estimations based on the model of aerogel in Sec. 12. II and 
the NMR data suggest that in the aerogel samples under investigation, the LIM 
length Llim at which the LROO is destroyed is of order of a micron. 

Though Llim is much bigger than the microscopic scales of aerogel strands, 
the orientational disorder may have pronounced effects. In particular, it may be 
the reason of existence of a small region of equilibrium A-phase in zero magnetic 
field 10 . For NMR experiments it is important that Llim is smaller than dipole 
length ^£> ~ 10 /im characterizing spin-orbit coupling between 1 and d in Eq.©. 
This leads to anomalously small values of the observed transverse NMR frequency 
shift and longitudinal NMR frequency— as we shall discuss in Sec. [3] 

A regular anisotropy may suppress the LIM effect, and this does take place 
when one applies a uniaxial deformation to the aerogel sample (Sec. l2.2l l. A linear 
deformation of about 1 % leads to the restoration of the uniform orientation of the 
1- vector in the sample 4 . As a result the value of the transverse NMR frequency 
shift is enhanced by an order of magnitude compared to that in the LIM state. If 
the deformation leads to orientation of 1 along the magnetic field H, the transverse 
NMR frequency shift becomes negative. These are the observations which support 
the interpretation of the A-phase like state in aerogel as the disordered LIM state. 

Open problems related to the disordered LIM state including the superfluid 
properties of the LIM state are discussed in Sec. [4] Superfluidity in LIM state could 
be rather unusual because, according to the Mermi-Ho relation in Eq.OJ between 
1 and \ s , the LIM texture generates superfluid vorticity. As a result, the disordered 
LIM state can be represented as a system of randomly distributed skyrmions - 
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vortices with continuous cores. The characteristic distance between the vortex- 
skyrmions and the size of their cores are determined by the LIM length Llim- In 
principle, the superfluid density of 3 He- A in aerogel may depend on the ability of 
the aerogel to pin the randomly distributed vortex-skyrmions. 



2 Theory of Larkin-Imry-Ma effect 

2. 1 Larkin-Imry-Ma state in the model of random cylinders 

We use a simple model of aerogel as a system of randomly oriented cylindrical 
strands (see Refs. n > 12 and Fig. [3] middle) with diameter 8 ~ 3 nm and length 
^ a ~ 20 nm, which is also the distance between strands. Since 8 is much smaller 
than the superfluid coherence length ^o, the theory of Rainer and Vuorio for the 
microscopic body immersed in 3 He-A 13 can be applied. According to 13 the 1- 
vector remains uniform outside the body (see Fig. [2] left), and the orientational 
energy acting from 1 on a cylinder with the direction of axis along n, is 

Ei = E a (t ■ n,) 2 , E a ~ ^4^5 > . (4) 

Here kp is Fermi momentum; the parameter E d is positive as follows from 13 . 

For the infinite system of cylinders, the average orienational effect of the ran- 
dom cylinders on the 1-vector is absent, if the 1-vector is kept uniform. However, 
there is a collective effect of many cylinders, which makes the 1-vector inhomo- 
geneous on large scales. Let us consider the box of size LxLxL, which contains 
a large but finite number N ~ L 3 3 ^> 1 of randomly oriented cylinders and still 
has a uniform orientation of the 1-vector. Due to fluctuations of orientational ener- 
gies of randomly oriented cylinders, the 1-vector will find the preferred orientation 
in the box, with the energy gain proportional to N x l 2 . The corresponding negative 
energy density is determined by the variance of the energy of random anisotropy: 

E ran r,-(T(E i -(E)) 2 \ L-*~-NV 1 Ejr*~-E&WirW. (5) 

The neighboring boxes prefer different orientations of 1, as a result the effect of 
orientation by the aerogel strands will be opposed by the gradient energy £ ?ra£ / ~ 
K{Vl) 2 ~ K/L 2 . The optimal size of a box with a uniform orientation of 1 within 
the box - the LIM length Llim at which LROO is destroyed - is found from 
competition between the energy of random anisotropy and the gradient energy: 

Eran + E grad - ~E a ^ 3/2 L~ 3 / 2 + KL~ 2 . (6) 
Minimization of Eq.© gives 

LuM ~ ~e\~ ~ ^t 2 " ' (7) 

Here we used Eq.© for Z? a ; and estimated the rigidity of 1 as K ~ (k 3 F /m) {A 2 /T 2 ). 
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Fig. 3 Uniaxial deformation of cylindrical sample of aerogel. Squeezing (stretching) of cylinder 
along its axis z leads to global anisotropy for 1 with easy axis along z (easy plane normal to z) 



Eq.© is valid for thin strings, 5 <C <5jo- This corresponds to collective pinning: 
1 is not pinned by individual strings, has a short-range orientational order and is 
perfectly uniform at the aerogel scale <!j a . The LROO is lost at scale Llim S> 
due to collective effect of many strings (Fig.^righi). With ^ a ~ ~ 20 nm and 
<S ^ 3 nm, the length Llim ~ 1 M m - Strong pinning would occur in case of thick 
strings, 8 > §a, when 1 is oriented by each string due to boundary conditions for 1 
at the surface of a string, and the LROO is destroyed at the scale of § a . 



2.2 LIM effect killed by regular anisotropy 

Squeezing of a cylindrical sample of aerogel along its axis (m, z ~ Al/l < 0, Fig. 
[3] /<?/?) makes z the easy axis for 1, while stretching (u zz ~ Al/l > 0, Fig. |3] rigfe) 
produces the easy plane anisotropy. In general, the deformation of aerogel 
induces regular anisotropy for the 1-vector with the density of orientational energy: 

E reg = (Ei - (E) ) \ = K (utk - ^Un Sg^j ljl k , fC - E^ 3 . (8) 

Comparing this energy with the LIM energy (in Eq.© at L = LlimX one finds that 
the ordered state with 1 oriented along the easy axis z becomes more energetically 
favorable than the LIM state when the squeezing is sufficiently large: 

I (^lim) ^ ^ 

In the weak pinning limit of 3 He- A in aerogel the LIM effect is so subtle, that 
the critical value of the deformation at which the first order phase transition from 
the LIM state to the uniform state occurs is rather small; Al/l - 1(T 3 - 10~ 2 . 
For large stretching the polar phase is predicted to appear in vicinity of r ca 14 . 
Experimentally, squeezing by about 1% is enough to get the uniform 1 texture 4 . 
On the other hand, recent experiments 15 in which compression has been applied to 
compensate intrinsic deformations, indicate that when the deformation becomes 
below 1%, the initial uniform 1 texture is transformed to the disordered state. 
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Fig. 4 left: Temperature dependence of positive frequency shift for two lines observed in NMR 
spectrum on 3 He-A in aerogel sample according to Ref. 3 (P = 29.3 bar, H = 224 G). right: Neg- 
ative frequency shift in the NMR spectrum observed in 3 He-A in deformed aerogel according to 
Ref.i (P = 29.3 bar, H = 290 G) 



3 Experiment 

3.1 Estimation of Larkin-Imry-Ma length from NMR 

Fig. [4] demonstrates results of the NMR experiments on 3 He-A in conventional 
aerogel sample 3 (left), and in a squeezed aerogel 4 (right). In the latter, the trans- 
verse NMR line has negative frequency shift corresponding to the uniform 1- vector 
oriented along the magnetic field H || z. 

In Ref. 3 two NMR lines have been observed in the transverse NMR spectrum 
in Fig. |4](/e/f)- The line with a larger frequency shift, denoted as /-line, can be 
removed after application of the 180° pulse while cooling through T ca (the su- 
perfluid transition temperature in aerogel). The salient feature of the c-line which 
survives after the 1 80° pulse is the anomalously small frequency shift. It is about 
30 times smaller than the magnitude of the negative frequency shift observed in 
Refs. 4 ' 5 in the aerogel sample with the uniform 1-vector (we take into account that 
the frequency shift is °c l/H with H = 224 G in Ref. 3 and H = 290 G in Ref.- 4 -). 
This large difference can be easily understood if one identifies the ostate (the state 
with a single c-line in the spectrum) as the LIM state. 

Properties of NMR on the LIM state depend on ratio between Llim and the 
dipole length ~ 10 Jim of spin-orbit coupling between 1 and d. The orbital vec- 
tor 1 will be locked with the spin-space vector d if Llim > £d an d unlocked from 
d if Llim < %D- The NMR results are in favor of the almost completely unlocked 
case, since in the limit Lum/^d *C 1, the 1-texture is fully disordered, (l 7 ;^ = 
1/3, and the frequency shift is zero 16 . The nonzero but relatively small magni- 
tude of the frequency shift compared to the negative frequency shift, A <D c _r me ~ 
zi conformed 1/30, can be ascribed to the (Lum/^d) 2 correction. This allows us to 
estimate Llim ~ 1 l^m, which is in a reasonable agreement with Eq.©. 
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Fig. 5 Tipping angle dependence of the frequency shift of transverse NMR in the c and / states 
of 3 He-A in aerogel according to Ref. 3 Experiments with compensated deformation demon- 
strated the ( 1 + cos /3 ) behavior similar to the /-line 15 



3.2 Tipping angle dependence of NMR 

Lines c and / in the spectrum of transverse NMR in 3 He-A in aerogel demonstrate 
different dependence of frequency shift on the tipping angle j3 in Fig. [5] This 
can be compared with the theoretical frequency shift in the limit of the almost 
completely disordered state, when Lum/^d *C 1: 

A( °(P) = & CO sj3+0(l + cos)3) . (10) 

4 ©max 

Here A(O max is the maximum possible frequency shift occurring in the uniform 
1-field, it coincides with the magnitude of the negative frequency shift when 1 || H; 
the parameter a <C 1 and id are nonzero due to the deviations from the full 
disorder caused by spin-orbit interaction—: 

i=i(l-3</?» , fl =i(l-2<sin 2 <£» , (11) 

and <P is the angle between the components lj_ and d^, which are transverse to 
magnetic field H || z. Parameters a and b are zero in the limit when disorder is 
strong compared to the spin-orbit interaction, i.e. when Lum/^d —> 0, and are 
proportional to (£lim/£d) 2 with pre-f actors depending on the type of disorder. 

The disorder can be described in the following phenomenological way. In the 
c-state the transverse component lj_ of the vector 1 is random, while the spin-space 
vector d is perpendicular to H and is regular. This gives a = 0. On the other hand 
the magnetic field H || z, via the vector d, produces a weak easy-plane anisotropy 
for the 1-vector. As a result, one has b > 0. The /-line represents regions with the 
/-state where the 1-vector is fully random, and thus b = 0, while d and lj_ are not 
fully independent, producing a > 0. 

However, it is not clear what is the microscopic background for such behav- 
ior of the parameters, and the detailed numerical simulations of the structure of 
different possible disordered states of 3 He-A in aerogel are required. Another in- 
terpretation of the fi -dependence of the two NMR lines in terms of the modified 
robust phase has been proposed by Fomin 17 . 
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4 Discussion: superfluid properties of LIM state 

In chiral anisotropic superfluids the LIM effect should influence the superfluid 
properties of the system. According to the Mermin-Ho relation in Eq.©, circu- 
lation of superfluid velocity along the closed path L is expressed in terms of the 
surface integral over the 1-field 




Here we normalized the circulation to the circulation quantum in 3 He, K= nh/m; 
the circulation number N is not necessarily integer in chiral superfluids. 

If we ignore for a moment that the 1-field is the part of the order parameter 
triad m, n and 1 in Eq.(fl]) and consider 1 as completely independent field, then the 
variance of N for sufficiently large contours is 

[2 

(At2( L ))^ , L»Llim. (13) 

LIM 

This corresponds to randomly distributed skyrmions - vortices with continuous 
cores of the size Llim, which is also the mean distance between the skyrmions. 
Such a state resembles the plasma of vortices above the Berezinskii-Kosterlitz- 
Thouless transition, where vortices destroy superfluidity. This suggests that in 
aerogel, while the local superfluid density p s on the scales below Llim ls of the 
same order as in 3 He-B, it could be reduced and even become zero at large scales. 

Experiments on 3 He- A in aerogel 1819 demonstrate that p s ^ though it is 
somewhat smaller as compared to p s in 3 He-B. This shows that in the consider- 
ation of the LIM effect, we cannot ignore the energy of superflow - the energy 
related to the gradient m • Vn, which is the superfluid velocity according to Eq.®. 
It was stated that the energy of superflow is not relevant for the LIM effect—. 
However, it should suppress the abundance of vortex-skyrmions. The pinning of 
skyrmions would also restore superfluidity. 

There are many open questions related to the LIM state. In particular, how rigid 
is the 1-texture at large length scales and how strongly are the vortex-skyrmions 
pinned by aerogel? The problem of rigidity in the LIM state is thirty years old, see 
the paper by Efetov and Larkin 21 and recent publications 22 . The idea of quasi- 
long-range order (QLRO) with a power-law decay of correlators discussed for a 
vortex lattice in superconductors 23 is applicable to 3 He-A in aerogel. Probably 
the QLRO does not occur in isotropic aerogel because of the non-Abelian SO(3) 
group of triad m, n, 1, but it should take place in a stretched aerogel where the easy 
plane anisotropy reduces SO(3) to its Abelian U(l) x U(l) subspace. Following 
Ref . 24 one may expect that the power law is determined by relations between all 
7 stiffness parameters K in the 3 He-A gradient energy 6 . The phase transition from 
LIM state with QLRO to the disordered LIM state or to the non-superfluid LIM 
state can be generated by increasing the density of skyrmions. The states c and / 
in NMR experiments 3 and states in hydrodynamic experiments 19 (see below) may 
differ by density of skyrmions and thus may have different superfluid properties. 

The external superflow influences the LIM state due to anisotropy of the su- 
perfluid density in Eq.©: since p s \\ < p s ±, the applied \ s produces the global 



9 




T/T c velocity mm/s 



Fig. 6 Anisotropic superfluid density of 3 He-A in aerogel according to Ref. 19 . Superfluid den- 
sity is history dependent and depends on the applied superfluid velocity. 



easy-axis anisotropy for 1. This suggests that the sufficiently large \ s destroys the 
LIM state, which is consistent with hydrodynamic experiment 19 demonstrating 
glassy behavior of 1-texture at low \ s and more regular behavior at large v^. If the 
initial state A in Fig. [6] is fully disordered but rigid, its superfluid density must be 
p^ = (1 /3)p s || + (2/3)p s ±. With increasing v, the state becomes more and more 

uniform. In the state H, which is reached at large v s , the vector 1 is oriented along 
y s , and one has pf = p s » < pf. Because of the hysteresis, the uniform texture 

persists when velocity is reduced, but now with 1 _L \ s as dictated by geometry of 
the sample; as a result in the subsequent state I one has p' = p s _i > pf. 

Another open problem is the role of extended objects, (i) Topological defects 
with hard cores (single-quantum and half-quantum vortices) may be pinned by 
aerogel. Trapping of singular vortices by aerogel has been observed in 3 He-B 5 . 
This suggests that c and / states may have different abundance of defects, (ii) The 
network of connected aerogel strands may have the long-range correlations (e.g. 
due to elastic deformations of aerogel). The interplay of the quenched long-range 
and short-range disorder may also lead to two LIM states with different scales^ 5 -. 

The LIM effect is rather subtle in 3 He-A, it is destroyed by weak anisotropy 
produced by deformation of aerogel and/or by flow. This suggests that maybe the 
adequate description of the LIM effect is in terms of statistics of random deforma- 
tions Ujk of aerogel with the energy term in Eq.®, instead of the discussed statis- 
tics of random orientations of strands. On the other hand, the global anisotropy 
produced by sufficiently large deformations allows us to study novel phenomena 
which cannot be observed in bulk 3 He- A. For example, orientation of 1 along the 
magnetic field H (not possible in a pure bulk 3 He-A) stabilizes: (i) Alice strings 
(half-quantum vortices) in rotating sample 26 ; and (ii) the phase-coherent preces- 
sion of magnetization 27 recently observed in a squeezed aerogel 28,5 . The coherent 
precession in 3 He-A is another realization of Bose-Einstein condensation (BEC) 
of magnons. The first example of magnon BEC found in 1984 in 3 He-B is known 
as the homogeneously precessing domain (HPD, see review s 29 ' 30 ), and histori- 
cally this was the first BEC which was experimentally stabilized. Magnon con- 
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densates in 3 He-A and condensates in 3 He-B (HPD, the state HPD2 observed in a 
stretched aerogel 31 and Q-balls 32 ) all have different spin-superfluid properties^. 

In conclusion, 3 He-A in aerogel is one of the most interesting objects for in- 
vestigation of the Larkin-Imry-Ma and related effects. 
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